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1.  Introduction 


The  mechanics  of  a  continuum  are  commonly  described  by  relating  two  sets  of  coordinates:  a  set 
of  reference  coordinates  that  serves  to  label  the  particles  in  an  arbitrary  (perhaps  initial  or  even 
fictitious)  configuration,  and  a  set  of  spatial  coordinates  that  fixes  locations  in  space  (that  is,  in 
the  laboratory)  (1 ,  2).  The  importance  of  the  reference  coordinates  in  continuum  mechanics  is 
that  they  serve  to  name  each  individual  particle  so  it  can  be  tracked  in  the  spatial  coordinates. 

This  is  the  sense  in  which  the  reference  coordinates  may  even  be  fictitious — since  they  merely 
serve  to  label  each  particle;  reference  coordinates  themselves  need  not  have  any  particular 
geometric  meaning  so  long  as  they  fulfill  this  requirement.  In  this  report,  we  seek  to  clarify  the 
relationship  between  the  two  standard  coordinate  sets  and  a  third  set  of  coordinates,  generally 
called  convective  coordinates  which  are  related  in  subtle  ways  to  both  the  spatial  and  reference 
coordinates.  Convective  coordinates  are  particularly  important  with  respect  to  nonmechanical 
physics  occurring  against  the  background  of  a  deforming  body,  and  are  almost  essential  within  the 
framework  of  the  theory  of  relativity,  as  that  theory  does  not  recognize  the  universality  of 
simultaneity. 

More  specifically,  in  this  report,  we  clarify  the  relationships  between  these  three  coordinate 
systems,  and  use  the  convective  coordinate  system  to  formulate  the  classical  theory  of 
electromagnetism  in  the  presence  of  material  deformation.  This  reformulation  is  useful  in 
numerical  work  because  the  physically  required  continuity  of  electromagnetic  field  descriptions 
depends  intimately  on  the  geometry  of  the  boundary,  which  is  generally  simple  to  describe  only  in 
reference  or  convective  coordinates.  By  reformulating  the  Maxwell  equations  in  convective 
coordinates,  we  also  demonstrate  that  the  standard  vector  formulation  of  them  is  not  covariant  in 
classical  physics  (3).  In  other  words,  we  demonstrate  that  a  proper  formulation  of  Maxwell’s 
equations,  implying  the  invariance  of  physical  law  for  all  observers,  cannot  be  achieved  within  the 
confines  of  classical  physics.  This  difficulty  is  a  primary  source  of  the  conflicting  formulation  of 
the  Maxwell  equations  in  continuum  mechanics  literature  (3-6). 


2.  Background  and  Curvilinear  Coordinates 


To  understand  the  subtle  differences  in  the  various  coordinate  systems  presented  in  this  work, 
definitions  need  to  be  made  clearly  and  with  a  modicum  of  rigor.  Moreover,  the  formulations 
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presented  in  the  following  concentrate  on  the  use  of  general  curvilinear  coordinates,  a  subject 
likely  unfamiliar  to  many  readers.  Therefore,  in  the  interest  of  precision  and  clarity,  we  describe 
curvilinear  tensor  theory  in  the  next  two  sections.  This  section  concentrates  on  the  algebra  of 
curvilinear  systems,  their  metric  structure,  and  the  vector  bases  used  to  define  field  quantities. 

The  next  section  describes  the  formulation  of  differential  operators  for  the  differentiation  of  field 
quantities  defined  in  such  spaces. 

2.1  The  Background  Cartesian  System 

Before  defining  the  physical  coordinate  systems  at  the  heart  of  this  report,  we  first  describe  an 
absolute,  mathematical  coordinate  system  disconnected  from  any  material  continuum  and  fixed 
firmly  as  a  background  description  of  space.  This  background  space  is  assumed  Cartesian  and 
fixed  for  all  time,  and  is  denoted  by  s.  The  space  s  is  constructed  given  a  point  O  (the  origin )  and 
three  orthonormal  vectors  ul5  u2,  and  u3.  An  arbitrary  point  P  has  coordinates  (&1,  x2,  x 3)  if  the 
vector  from  0  to  P  is  given  by  the  position  vector 

x  =  XlVii.  (1) 

Superscripts  are  used  here  for  indexing  for  reasons  that  will  become  clear  later.  In  any  case,  here 
we  use  the  Einstein  convention  that  an  index  repeated  in  a  superscript  and  a  subscript  is  to  be 
summed  over  its  range. 

In  addition  to  serving  as  the  basis  for  geometrical  description,  s  can  be  used  to  describe  physical 
scalar  and  vector  fields.  For  instance,  if  we  have  a  (static)  pressure  associated  with  each  point  in 
space,  we  may  write 

P  =  p(xl),  (2) 

where  the  xl  refer  to  points  of  space  as  defined  by  equation  1 .  Similarly,  the  electric  field  at  the 
point  with  coordinates  xl  can  be  written  in  the  form 

e  =  e?(V)u  u  (3) 

where  now  both  the  coordinates  and  the  basis  vectors  of  underlying  background  description  are 
employed. 

2.2  Curvilinear  Coordinates  and  Metrics 

Our  initial  description  is  based  on  Cartesian  coordinates  because  such  a  description  allows  us  to 
define  coordinates  through  equation  1  and  thus  directly  connect  coordinates  with  position  vectors. 
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On  the  other  hand,  any  set  of  three  numbers  (x1,  x2,  x3)  suffices  to  locate  points  in  space  if  every 
such  triplet  is  mapped  uniquely  to  a  point  in  space;  that  is,  if  the  mapping 

^  =  Me  {1,2,3}  (4) 

has  the  properties  that  each  coordinate  set  xl  is  the  image  of  a  exactly  one  coordinate  set  x\ 

(From  this  point  on,  the  range  on  indices,  which  for  this  work  is  always  {1,  2,  3},  is  assumed  and 
suppressed.)  That  is,  we  are  assuming  the  mapping  is  bijective,  i.e.,  there  is  a  unique  set  of 
functions  such  that 

^  =  e(x*).  (5) 

These  new  x1  coordinates  function  as  well  for  point  location  as  did  the  old  x'  coordinates  because 
of  this  one-to-one  mapping:  Given  the  x\  the  xl  are  determined  uniquely,  and  associated  with  the 
point  described  by  equation  1.  We  therefore  refer  to  this  new  curvilinear  coordinate  systems  as  s, 
and  label  all  vectors  in  5  just  like  their  counterparts  in  s,  but  with  a  tilde.  The  only  weakness  of 
using  5  to  describe  space  is  that  nothing  as  simple  as  equation  1  relates  the  point  location  to  the 
coordinates. 


Given  this  failure  of  equation  1  and  the  unclear  relationship  between  the  x1  and  the  u,,  it  makes 
sense  to  define  new  basis  vectors  to  describe  vector  fields  in  s  (if  not  position  vectors).  Inspired 
by  the  original  Cartesian  basis  vectors,  we  might  define  our  new  basis  to  be  in  the  direction  of 
increase  of  a  single  coordinate,  computed  holding  the  other  coordinates  constant.  Specifically, 
define 


(6) 


These  vectors  are  not  necessarily  of  unit  length  or  orthogonal,  but  they  do  point  in  the  direction 
desired.  The  basis  is  illustrated  at  a  fixed  point  in  figure  la.  These  vectors,  and  any  other 
tensorial  quantities  indexed  by  a  subscript  and  therefore  subject  to  variable  changes  in  the  manner 
of  equation  6,  are  called  covariant. 


The  Cartesian  basis  vectors  are  not  only  well  suited  to  their  coordinate  system  because  they  point 
in  the  direction  of  increase  of  a  coordinate  holding  all  other  coordinates  constant;  they  are  also 
orthogonal  to  the  constant  coordinate  value  surfaces  constructed  by  holding  a  single  variable 
constant  and  letting  all  other  coordinates  vary.  The  u*  do  not  have  this  property.  In  view  of  this, 
we  can  define  a  new  set  of  basis  vectors 


u 


l 


dx  , 

— — rlT 

dxl 


(7) 
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Figure  1.  (a)  Covariant  basis  vectors  u*,  and  (b)  contravariant  basis  vectors  u*  at  the  point  (b,  b).  The  grid 
is  drawn  assuming  a  <  b  <  c  <  d  . 


where  for  notational  convenience  we  define 


u 


i 


(8) 


and  as  usual  S'3  (with  indices  as  superscripts  or  subscripts  or  both)  is  the  Kronecker  delta.  These 
vectors  are  orthogonal  to  the  constant  coordinate  value  surfaces  (i.e.,  fi'  is  orthogonal  to  the 
x 1  =  constant  surface)  since  they  are  the  gradients  of  the  curvilinear  coordinate  values  with 
respect  to  the  underlying  Cartesian  system.  Indeed,  the  biorthogonality  of  the  covariant  and 
contravariant  components  follows  from  the  chain  rule: 

_  fdx1  A  /  dx3  \  „  dxl  dx3  dx3  dx3  - 

U  Uj  \  dx1  J  \  dx3  3 )  3  dx1  dx3  did  dx3  3 

They  are  illustrated  in  figure  lb,  and  are  called  contravariant  because  they  change  in  a  manner 
opposite  that  of  the  covariant  basis  as  expressed  by  equation  7. 


Finally,  we  note  that  we  can  define  higher-order  tensors  that  obey  similar  rules,  and  that  they  can 
be  covariant  or  contravariant  in  any  of  their  indices.  In  this  work,  no  tensors  of  higher  than 
second  order  appears.  As  an  example,  a  twice-contravariant  set  of  second-order  tensor 
coefficients  ad3  transforms  according  to  the  formula 


dx  dx3  - 
a3  =  -tt-tt-cl3 
dx 1  dx3 


(10) 
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precisely  because  the  vector  outer  product  Ujir,  is  twice-covariant: 


dx1  dxj 


U^  =  ^^UiU- 


(ID 


Indeed,  these  two  equations  taken  together  ensure  the  product  of  such  quantities  is  invariant,  since 


a*Ju*uj  = 


dx1  dxJ 
dxm  dxm ' 


=  5lmS3namn  UjUj 


=  alJUiUj. 


(  dxl  dxj  \ 
\  dx1  dxdU,U> ) 


(12) 


This  invariance  under  the  change  of  coordinates  is  a  hallmark  of  a  correctly  formulated  physical 
theory  since  it  ensures  that  the  meaning  of  physical  quantities  is  independent  of  their 
mathematical  description.  Similarly,  a  mixed  second-order  tensor  a*  changes  coordinates  in  the 
by  now  expected  fashion 


dx 1  dxj  j 
■1  dx1  d'Jd 


(13) 


because  the  product  u,iiJ  changes  in  the  opposite  fashion. 


Because  of  these  notational  observations,  it  is  easy  to  compute  the  inner  product  of  two  vectors 
expanded  with  respect  to  the  two  different  bases.  If  we  write 


a  =  a*Uj  =  aju*, 


(14) 


and 


b  =  6*Uj  =  bi  u*, 


then 


a  •  b  =  dlbi  =  a^b1. 


(15) 

(16) 


(As  an  aside,  note  that  the  boldface  vector  notation  does  not  distinguish  between  covariant  and 
contravariant  coordinate  systems.  This  is  because  vectors  themselves,  as  independent,  physically 
meaningful  entities,  are  by  definition  independent  of  coordinate  system.  Basis  vectors,  which 
seem  to  be  an  exception  to  this  rule,  are  in  fact  not  exceptions:  For  them,  the  use  (or  nonuse)  of  a 
tilde  indicates  for  which  space  they  form  a  covariant  or  contravariant  basis,  rather  than  having 
anything  to  do  with  their  expansion  in  any  system.) 


To  ease  the  computation  of  scalar  products  between  pairs  of  covariant  or  contravariant  vectors,  a 
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metric  tensor  is  introduced.  Defining  the  metric  tensor 


9i~]  =  u*  ■  Uj,  (17) 

the  inner  product  of  a  and  b  can  be  computed  from  contravariant  components  alone: 

a-b  =  cjijtfV .  (18) 


In  a  similar  manner,  the  twice  contravariant  metric  tensor  is  defined: 

gij  =  u*  •  uJ.  (19) 

These  two  tensors  are  trivially  symmetric  and  are  inverses  of  one  another  by  the  chain  rule: 


(20) 


The  introduction  of  the  metric  tensor  creates  the  possibility  of  working  entirely  with  vector 
coefficients,  and  assuming  and  suppressing  the  basis  vectors  that  are  always  defined  by 
equations  6  and  7.  (Indeed,  some  books  (7)  never  even  mention  these  vectors,  assuming  them 
superfluous.)  For  this  approach  to  be  useful,  the  metric  tensor  should  be  computable  from  the 
functions  £  and  £,  and  indeed  it  is.  The  standard  covariant  metric  tensor  may  be  computed  from 


Hij  Uj  •  Uj 


5, 


dx * dxj 
lj  dx 1  dx- ' 


(21) 


recalling  in  the  final  step  that  s  is  orthonormal  so  u,  •  uv  =  .  This  computation  has  the  added 

benefit  of  demonstrating  that  is  a  tensor,  since  the  metric  tensor  in  s  is  trivially  given  by 

=  Sij.  In  precisely  the  same  way,  the  inverse  of  the  metric  tensor  g13  can  be  directly  computed 
according  to  the  formula 


au  Sj'  —  — 

dx{  dxj  ’ 


(22) 


and  is  a  twice  contravariant  tensor. 


Finally,  computing  the  dot  product  of  the  vector  a  defined  above  with  each  of  the  basis  vectors  in 
turn  demonstrates  that  the  covariant  and  contravariant  components  of  a  vector  are  related  through 
the  metric  tensor: 

cii  gjj(d .  (23) 
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and 


(24) 


a1  =  gtJaj. 


These  operations  are  known  as  the  lowering  and  raising  of  an  index,  respectively. 

2.3  Determinants  and  Cross  Products 

In  addition  to  the  computation  of  dot  products  (and  thereby  lengths),  work  in  physics  requires  the 
computation  of  cross  products  and  determinants.  These  computations  depend  intimately  upon 
the  definition  of  the  Levi-Civita  system 

0  if  i  —  j,  j  —  k,  or  i  —  k, 

etjk  —  \  1  if  ij k  is  an  even  permutation  of  123,  (25) 

-1  if  ij  k  is  an  odd  permutation  of  123. 

This  definition  shows  clearly  that  all  completely  antisymmetric  third-order  systems  are  scalar 
multiples  of  the  Levi-Civita  system,  a  fact  we  will  need  later.  We  can  also  define  the 
contravariant  Levi-Civita  system  in  the  obvious  way,  specifically 


eijk  =  6irSjsSkterst 


(26) 


Determinants  and  Jacobians  can  be  easily  computed  in  terms  of  this  tensor.  The  Jacobian  of  the 
mapping  from  the  xl  to  the  xl  is  given  by  the  formula 


det 


dxl 

dx1 


^-ijk 


dx1  dx^  dxk 
dx 1  dx 2  dx 3 


(27) 


This  observation  can  be  used  to  clarify  the  tensorial  nature  of  the  Levi-Civita  system. 
Specifically,  we  can  consider  the  system 


^ ink  ~ 


dx1  dx^  dxk 


rjk  Qtf  d&  dxk' 

This  system  is  completely  antisymmetric  because,  for  instance, 

dx 1  dxi  dxk  dx^  dx 1  dxk 

xk  =  eiikd&dtfmt  =  -6jikdtfd&d5*  =  ifk' 

Therefore,  because  all  antisymmetric  systems  of  third  order  are  uniquely  determined  up  to  a 


(28) 


(29) 
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multiplicative  constant,  we  find  that 


dxl  dxj  dxk  (  dx1  \ 

€zjkdxl  dxi  dxk  C6^  \dxl) 


(30) 


and,  by  the  same  token, 


iikdxld^dxk  =  ^  ^ 

dxi  dxi  dxk  \  dxi ) 


(31) 


These  two  equations  imply  that  the  Levi-Civita  system  is  not  a  tensor  in  the  usual  sense,  but  in  a 
new  sense  in  which  the  usual  tensorial  transformation  is  accompanied  by  multiplication  with  a 
power  of  the  Jacobian  determinant.  Such  tensors  are  called  relative  tensors ,  or,  if  the  power  of 
the  determinant  is  ±1,  tensor  densities  (7).  (Note  that  no  tilde  is  ever  put  on  e  since  it  has  the 
same  values  in  all  systems  and  so  needs  no  such  distinction.)  In  this  sense,  we  may  refer  to  it  as 
the  “Levi-Civita  tensor.” 


Because  the  determinant  of  the  variable  transformation  tensor  occurs  so  frequently  in  the 
following,  we  can  simplify  the  remaining  exposition  with  better  notation.  If  we  let  the 
determinant  of  a  twice-covariant  tensor  be  represented  by  its  unsubscripted  symbol,  we  can  write 


g  =  det  (<y 


(32) 


or,  more  simply, 

'/F=det(S?)'  <33) 

Note  again  that  this  definition  pertains  only  to  the  determinants  of  twice-covariant  tensors,  and  in 
particular  to  the  twice-covariant  metric  tensor. 


This  clean  expression  for  the  Jacobian  can  be  used  to  simplify  the  expression  of  cross  products  in 
the  curvilinear  system.  If  the  original  Cartesian  space  is  right-handed,  as  we  shall  always 
assume,  the  cross  product  c  =  c*Uj  =  qiT  of  the  vectors  a  and  b  of  equations  14  and  15  has 
contravariant  components  given  by 

c*  =  eijkajbk  (34) 

and  covariant  components  given  by 

c%  =  <  (35) 


8 


Applying  all  the  change  of  basis  rules  described  above  to  these  formulas,  we  find  that  the 
coefficients  of  the  cross  product  are  given  in  the  curvilinear  system  by 


and 


c*  =  —f&afa,  (36) 

Ci  =  y/ge^apy*.  (37) 


3.  Differential  Operators  in  Curvilinear  Spaces 


Having  completed  the  description  of  the  algebra  underlying  the  description  of  fields  in  curvilinear 
coordinate  systems  in  the  last  section,  we  turn  to  the  differentiation  of  such  fields  in  this  section. 
The  main  idea  is  that  tensors  should  be  differentiated  in  a  physically  meaningful  way;  that  is,  that 
the  results  of  such  differentiations  should  be  tensors.  A  so-called  covariant  derivative  ensures  the 
physicality  of  results.  It  is  defined  immediately  below. 

3.1  The  Covariant  Derivative 


The  last  section  described  the  algebra  involved  in  changing  coordinates  between  a  fixed  Cartesian 
system  s  and  a  curvilinear  system  5.  In  this  section,  we  explore  the  analytical  properties  of 
coordinate  system  changes.  In  particular,  given  a  vector  a,  we  wish  to  determine  how  its 
derivatives  with  respect  to  space  can  be  computed  such  that  their  meaning  is  not  bound  to  any 
particular  coordinate  system. 

To  this  end,  consider  the  change  undergone  by  a  vector  a  over  a  differential  distance  dx.  In  the  5 
system,  the  change  in  a  can  be  expressed  in  terms  of  a  dot  product  as 

da  =  —  dx\  (38) 

dx* 


To  actually  apply  this  formula,  a  more  specific  formula  for  the  components  of  the  partial 
derivative  indicated  above  is  needed.  To  this  end,  we  define  the  covariant  derivatives  of  the 


components  as 


Da  j  <9a  _ 

- A  =  •  U;, 

Dad  dx 1  3 


(39) 


and 


D  a?  .  <9a  ~ 

— —  =  -  •  U. 

Dad  dx 1 


(40) 
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The  word  “covariant”  used  here  in  concert  with  the  derivative  definition  is  not  meant  to  indicate 
the  nature  of  an  index  (though  covariant  differentiation  does  result  in  the  addition  of  a  covariant 
index  to  a  tensor).  Rather,  the  word  implies  that  the  result  is  tensorial;  that  is,  it  is  covariant  in 
that  it  changes  form  appropriately  under  a  change  of  coordinates. 


To  find  practical  formulas  for  the  covariant  derivative  components,  we  begin  with  equation  39  and 
apply  the  product  rule  of  differentiation  to  find 


aa  asj.j  .  9u‘ 

a¥  ~  ~d¥  Wu  +  <41) 

Computing  the  inner  product  of  this  equation  with  the  basis  vector  uj,  and  defining  the  Christoffel 
symbol  of  the  second  kind  by 


k 


duk  r  (9u,- 

U;  •  — —  =  rC  •  -—dr, 

J  dxl  dxl 


(42) 


we  find  that  (7,  8) 


Doj  ddj 
Dxl  dxl 


k 


(43) 


Other  authors  use  the  notation 

1  ■  ! . }  <44) 
for  the  Christoffel  symbols,  but  this  notation  is  not  employed  here.  In  particular,  the  alternative 
notation  makes  the  Christoffel  symbols  look  like  third-order  tensors,  which  they  are  not. 


The  covariant  derivative  of  the  contravariant  components  of  a  vector  can  similarly  be  shown  to  be 
(7,8) 


Deri  _  d  aJ  I  j  Ijj 
D:r*  dxl  |  k  1  J  ° 

The  definition  can  even  be  extended  to  covariant,  contravariant,  and  mixed  tensors  of  second 


(45) 
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order  yielding  (7,  8) 


Dajj  da,j 


Jdxk 
D  a13 


dxk 

da13 


~  1 

1  1 

f  ~  ) 

m 

f  a«Xi  —  \ 

1  rn  1 

ik  J 

l  ~k~3  1 

~mj 


Dxk  dxk 


D  a% 


d~a j 


+ 


in  k 
i 

\  a'f 
ink  1 


Dxk  dxk 

as  the  appropriate  formulas  for  the  coefficient  derivatives. 

3.2  The  Christoffel  Symbols 


k 

in  j 

in 

kj 


Q'rrv 


(46) 

(47) 

(48) 


The  Christoffel  symbols  used  in  the  definition  of  the  covariant  derivative  and  defined  in 
equation  42  obey  some  important  identities  that  find  uses  in  this  work.  These  identities  are  for 
the  most  part  proven  by  complicated  (but  purely  formal)  manipulations,  and  their  proofs  can  be 
found  in  almost  any  work  on  tensor  analysis,  such  as  references  7  or  8. 

First,  it  can  be  demonstrated  that  the  symbols  are  symmetric  in  their  lower  indices;  that  is,  that 


so  that  in  particular 


(49) 


(50) 


This  symmetry  can  be  used  to  derive  their  expansion  in  terms  of  the  metric  tensor,  which  after 
some  tedious  manipulations  can  be  shown  to  be  (7,  8) 


i 

jk 


t  dtirnk 

dxk 


(51) 


Two  more  facts  about  Christoffel  symbols  are  important  in  the  remainder  of  the  work.  The  first  is 
Ricci’s  lemma,  which  states  that  the  covariant  derivative  of  the  metric  tensor  vanishes.  The  proof 
is  simple:  In  s,  gl3  =  d,3 ,  so  its  covariant  derivative  vanishes  trivially.  Since  the  metric  tensor  is 
itself  a  tensor,  its  covariant  derivative  is  also  a  tensor.  Therefore,  its  vanishing  in  one  system 
proves  its  vanishing  in  all  systems.  Ricci’s  lemma  implies  that  the  metric  tensor  can  be  treated  as 
a  constant  for  the  purposes  of  covariant  differentiation. 
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Finally,  intricate  (but  ultimately  straightforward)  manipulations  of  the  Christoffel  symbol 
definition  show  that  (7,  8) 


i  1  =  1  dyfg 

ij  I  VF  d&  ■ 


(52) 


3.3  Differential  Operators 

With  all  of  this  detail  in  hand,  the  standard  differential  operators  of  multivariate  calculus  can  now 
be  introduced.  In  s,  the  gradient  of  a  scalar  (f>(xl)  has  components  equal  to  the  partial  derivatives 
with  respect  to  the  underlying  coordinate  dimensions: 

<53) 

Note  that  the  components  of  the  gradient  are  covariant,  as  they  are  associated  with  the 
contravariant  basis  vectors  ensuring  the  consistency  of  the  Einstein  summation  convention. 
Assuming  that  0  has  a  physical  meaning  independent  of  the  coordinate  system  in  which  it  is 
defined  (that  is,  assuming  0  a  true  physical  scalar),  its  value  in  s  must  obey 

0(57)  =  4>(xl)  .  (54) 


To  find  the  gradient  value  in  s,  we  invoke  the  formula  for  contravariant  change  of  basis, 


<90 

dxi 


dcj)  f  dxl 
dxi  \  dx 1 U  )  ’ 


(55) 


so  that,  by  the  chain  rule,  the  gradient  is  given  by 


V0 


(56) 


The  divergence  of  a  vector  a  is  given  by  the  formula 


Va 


Da* 

Da;* 


(57) 


in  s.  Covariant  differentiation  is  used  here  to  ensure  that  the  result  is  a  true  scalar.  (It  is  trivially 
identical  to  partial  differentiation  here  anyway  given  the  Cartesian  assumption  on  5.)  Now, 
because  covariant  derivatives  change  coordinates  according  to  the  laws  of  tensor  algebra,  we  can 
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immediately  write 


V  •  a  = 


Using  equations  45  and  52,  we  can  expand  this  expression  to  find 


Da*  da1 
Da;*  dx1 


=^_+  1  dy/F, 

dx1  yfg'  dxi 


or  finally,  by  separating  a  factor  of  g  2  and  recognizing  the  derivative  of  a  product, 


v ' a  =  WIF 


Lastly,  there  is  the  curl.  In  s,  the  curl  is  defined  as 

V  x  a  =  e^k — -Uj. 

Dxi 

This  formula  can  be  brought  into  s  using  equation  31  and  the  usual  covariant  variable  change 
formulas  to  find 


Dxi1 


1  dxl  dxi  dxk 
yfg  dx ™  dxh  dxP  ' 


J_erhhpSmSnSP^k^ 

yj g  J  kDxi 

1  ifk^ak~ 

VF  - 


dxi  Qxk  Dak  \  /  dx'  _ 
dxi  ()xk  Dxi  )  V  dxi  U’ 


In  view  of  equation  50,  however, 


cw k^ak  _  cwk 


so  that  the  curl  may  finally  be  written  without  the  covariant  derivative  notation  as 


V  x  a  = 

yjg  dxi 
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4.  The  Coordinate  Systems  of  Continuum  Mechanics 


4.1  Reference  Coordinates 

To  begin  applying  the  mathematics  of  the  foregoing  sections  to  continuum  mechanics,  coordinate 
systems  must  be  introduced  to  describe  continua  and  their  motions  and  deformations.  The  first 
set  of  coordinates  introduced  in  this  regard  is  the  reference  coordinate  set ,  which  serves  to  name 
the  particles.  Indeed,  in  more  general,  mathematical  renderings  of  continuum  mechanics,  the 
reference  “coordinates”  hardly  even  need  to  be  coordinates  or  even  numerical,  so  long  as  each 
particle  is  identified  as  a  member  of  a  set  (9).  We  therefore  assume  that  a  material  body  is 
composed  of  particles,  which,  at  some  time,  could  have  been  identified  by  their  location  X1  in  a 
Cartesian  coordinate  system  (5.  By  the  seeming  circumlocution  “at  some  time,  could  have  been,” 
we  mean  to  imply  that  while  the  particles  may  never  have  been  in  the  reference  configuration  X1, 
the  mapping  from  that  configuration  to  any  future  configuration  xl  is 

•  Injective,  so  that  each  particle  in  the  current  configuration  corresponds  to  only  one 
reference  particle, 

•  Surjective,  so  that  each  particle  in  the  current  configuration  has  a  name,  and 

•  Characterized  by  a  mapping  with  a  positive  Jacobian  determinant,  so  that  the  chirality  of 
the  mapping  is  preserved  throughout  space  and  time. 

In  most  work  in  elasticity,  the  X1  are  imagined  to  represent  the  rest  state  of  the  matter,  devoid  of 
elastic  potential  energy.  In  other  work,  the  reference  coordinates  may  be  the  initial  state  of  the 
continuum,  or  they  may  even  be  ignored  altogether,  except  to  assume  their  existence.  (This  last 
approach  is  often  seen  in  fluid  mechanics.) 

Finally,  in  this  work,  we  always  assume  that  ©  is  Cartesian,  as  mentioned  above  but  this  need  not 
be  the  case  in  more  general  applications.  For  some  problems  involving  lower  dimensional 
formulations  of  the  theory  or  special  symmetry,  it  might  be  useful  to  introduce  curvilinear 
material  coordinates  through  the  mapping  X1  =  E^X1).  By  the  assumptions  listed  above  and 
the  discussion  of  the  previous  two  sections,  these  new  coordinates  would  name  the  particles  as 
well  as  the  original  coordinates,  and  could  therefore  be  used  as  the  basis  of  a  theory  with  no 
essential  complications. 
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4.2  Spatial  Coordinates 


The  current  state  of  the  particles  of  the  body  named  by  the  reference  coordinates  is  related  by  the 
spatial  coordinates  xl.  The  location  of  a  given  particle  X1  at  a  given  time  t  is  completely 
specified  by  the  motion  of  the  particle 


^  =  e(X7,f).  (65) 

This  relation  completely  specifies  particle  kinematics  and  is  the  basis  for  the  theory  of  continuum 
mechanics. 

At  any  fixed  time  t  =  t(h  this  relationship  must  be  bijective  with  a  positive  Jacobian  determinant 
as  mentioned  in  the  previous  section.  Because  of  this,  we  may  also  write 

X1  =  EI(xi,t)  .  (66) 

This  equation  should  not  be  interpreted  as  meaning  that  the  X 1  depend  on  time;  they  don’t.  For 
an  arbitrary  fixed  value  of  time  t  =  t0,  what  equations  65  and  66  mean  is  that  £  and  S  are  inverse 
functions,  i.e.,  that 

a;1  =  e(X7,f0)  =  e{^{x\t0)  ,t0)  ,  (67) 

or,  equivalently, 

x1  =  =  S7(e(X7,f0)  ,t0)  ,  (68) 

for  all  t0. 

Because  the  particles  move,  the  relationship  between  their  current  location  and  their  reference 
designation  does  evolve,  and  that  is  captured  by  equation  66.  When  combined  with  the 
curvilinear  coordinate  information  provided  in  the  previous  sections,  this  observation  can  be  used 
as  a  basis  for  the  definition  of  the  coordinate  system  at  the  center  of  the  work,  which  is  provided 
in  the  next  subsection. 

4.3  Convective  Coordinates 

The  description  of  boundaries,  inclusions,  and  other  geometrically  relevant  physical  features  of  a 
body  are  often  more  simply  described  in  the  reference  system  ©  than  in  the  spatial  system  s. 

After  all,  the  reference  system  is  merely  a  matter  of  particle  nomenclature,  whereas  the  particle 
trajectories  in  s  are  governed  by  the  applied  excitations  and  the  laws  of  physics.  This  can  be 
illustrated  through  the  example  of  an  elastic  bar  in  the  shape  of  a  rectangular  prism  whose  axes 
are  aligned  with  the  coordinate  axes,  subjected  to  a  large  force.  Assuming  that  the  reference 
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coordinates  are  taken  to  be  the  initial  coordinates  of  the  bar  at  rest  at  t  =  0,  this  situation  is 
depicted  in  figure  2a.  In  &,  the  boundaries  of  the  bar  correspond  to  equations  of  the  form 
X 1  =  K1,  where  K1  is  some  fixed  constant  for  each  boundary  in  each  dimension.  Of  course,  at 
t  —  0,  the  boundaries  in  s  are  just  as  simple  to  describe,  since  at  this  instant  they  coincide. 
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Figure  2.  (a)  The  initial  state  of  the  bar  labeled  in  the  spatial  (dashed,  s)  and 
reference  (solid,  ©)  coordinate  systems,  (b)  the  deformed  bar  la¬ 
beled  in  spatial  coordinates  (5),  and  (c)  the  deformed  bar  labeled 
in  convective  coordinates  (5).  The  reference  coordinates  never 
change,  remaining  associated  with  (a)  even  as  the  deformation  oc¬ 
curs  . 


At  a  later  time,  under  the  action  of  outside  forces,  the  bar  will  have  deformed.  Suppose  it  takes 
the  shape  depicted  in  figure  2b.  The  spatial  coordinates  are  depicted  in  this  picture,  and  the 
boundaries  (obviously)  do  not  align  with  them.  By  definition,  the  reference  coordinates  never 
change,  and  are  still  depicted  in  figure  2a. 

To  arrive  at  a  coordinate  system  in  which  features  simply  described  in  reference  coordinates 
retain  their  simple  description  as  the  body  evolves  and  deforms,  we  turn  to  a  set  of  spatial 
curvilinear  coordinates  commonly  called  convective  coordinates.  Since  they  are  curvilinear 
spatial  coordinates,  we  denote  the  convective  variables  by  xl  and  the  convective  system  by  s.  The 
mapping  that  gives  rise  to  these  coordinates  can  be  written,  at  each  instant  t  as 

xl  =  €  {x\  t )  =  5'jE1  (x\  t )  ,  (69) 
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where  the  function  S  is  defined  in  equation  66.  (Note  that  these  coordinates  evolve  as  a  function 
of  the  “parameter”  t .)  The  virtue  of  this  choice  is  that  at  any  arbitrary  time  t0,  the  particle  with 
material  coordinates  X1  in  ©  has  coordinates  x1  =  X1  in  s.  After  all,  by  equation  68, 

X*  =  to)  =  SjE^eiX1,  to)  ,  to)  =  fix1.  (70) 

In  particular,  the  algebraic  description  of  interesting  object  features  never  changes.  This  is 
depicted  in  figure  2c,  which  labels  the  deformed  bar  of  figure  2b  with  the  convective  coordinates 
of  system  s. 

The  system  $  does  both  of  the  following: 

•  Describes  the  system  as  it  evolves,  even  if  the  “continuum”  depicted  breaks  apart  during  its 
evolution,  and 

•  Tracks  the  trajectory  of  fixed  particles  during  that  evolution. 

This  can  be  very  useful  for  the  discussion  of  nonmechanical  physics  unfolding  against  the 
background  of  the  continuum.  In  particular,  because  simultaneity  is  not  absolute  in  the  theory  of 
relativity  (10,  11),  convective  descriptions  can  become  essential.  In  the  next  section,  operations 
necessary  for  the  formulation  of  continuum  kinematics  in  convective  coordinates  are  described. 


5.  The  Operations  of  Continuum  Mechanics  in  Convective  Coordinates 


5.1  Temporal  Differentiation 

Temporal  differentiation  in  s  is  not  very  different  than  in  &  because  at  all  times  x'  =  X\ 
Therefore,  many  of  the  formulas  presented  here  have  direct  reference  counterparts  discussed  at 
length  in  standard  texts  on  continuum  mechanics  such  as  references  1,  2,  and  9. 


Perhaps  the  most  important  temporal  derivative  definition  is  also  the  most  basic:  the  particle 
velocity,  which  is  computed  by  holding  the  material  coordinates  (the  X1)  of  the  particle  fixed,  and 
differentiating  the  spatial  coordinates  with  respect  to  time: 


v 


i  _i 


dad 

d  t 


(71) 


The  use  of  the  “total  derivative”  notation  d/d t  here  is  meant  to  indicate  that  the  X1  remain  fixed 
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but  the  xl  vary  during  the  differentiation.  The  resulting  vector  if  is  the  velocity  of  a  fixed  particle 
through  space. 


Another  type  of  velocity  that  can  be  measured  is  the  speed  with  which  particles  flow  past  a  given 
point  in  space.  To  compute  this,  we  need  to  fix  the  point  in  space  (the  x1),  and  differentiate  the  x1 
with  respect  to  time.  We  can  thus  define 


ft 


dxl 

~dt' 


(72) 


where  the  “partial  derivative”  notation  d/d t  indicates  that  the  xl  remain  fixed  and  the  X1  are 
allowed  to  vary  as  particles  pass.  The  definitions  of  the  temporal  differentiation  operators  above 
should  make  it  clear  that 


(73) 


and 


(74) 


Finally,  functions  are  often  presented  in  terms  of  their  dependence  on  the  spatial  variables,  but 
need  to  be  differentiated  holding  the  particle  constant.  Given  a  function  /  =  f(xl,t),  its  material 
derivative  is  computed  as 

^l  =  ?L  +  ^L^L  =  ?L  +  vi?L  (15) 

dt  dt  dxi  dt  dt  dxi 


5.2  Deformations,  Transformations,  and  Metrics 

Consider  two  points  in  s  separated  by  an  infinitesimal  displacement  dx\  This  displacement  is 
related  to  an  infinitesimal  displacement  dX !  in  ©  by  the  equation 

dxl  =  g^jdX1,  (76) 

and  to  an  infinitesimal  variable  transformation  dxl  in  s  by 

dxl  =  —  dx\  (77) 

dxl 

Equations  of  this  form,  of  course,  hold  for  any  change  of  variables;  after  all,  this  is  just  a 
straightforward  application  of  the  chain  rule  of  calculus.  Moreover,  because  X1  =  x1  at  all 
times,  the  set  of  numbers  invoked  on  the  right-hand  sides  of  equations  76  and  77  are  the  same. 
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Nonetheless,  their  meanings  are  different  because  they  refer  to  spaces  with  different  metric 
structures.  The  squared  magnitude  of  dA"7  is  given  by 

dS2  =  6udXrdXJ,  (78) 

whereas  the  squared  magnitude  of  dxl  is 


d  s2  =  gij  dxldxP.  (79) 

The  reason  for  this  “discrepancy”  is  that  two  different  lengths  are  involved:  dS2  is  the  square 
length  of  the  reference  segment  that  corresponds  to  dx  in  the  spatial  domain,  whereas  ds2  =  ds2 
because  they  are  the  squared  length  of  the  same  segment  dx  in  different  coordinate  systems.  Of 
course,  by  equation  76  we  may  write 

ds2  =  C/jdX7dXJ  (80) 


where  the  Green  deformation  tensor  (7)  or  right  Cauchy-Green  tensor  (72)  is  defined  by 

dx 1  dxj 
13  =  bt2dXT  dX  J 


(81) 


Recalling  that  x1  =  X1  at  all  times,  and  recalling  the  definition  of  the  metric  tensor  in  terms  of 
the  variable  transformations  given  in  equation  21,  we  find  that 


9i-i  =  HSjCij, 


(82) 


that  is,  that  the  right  Cauchy-Green  tensor  and  the  convective  metric  tensor  are  numerically 
indistinguishable. 


From  these  observations,  the  computation  of  (oriented)  areas  and  lengths  using  reference, 
convective,  and  spatial  coordinates  is  simple.  A  differential  oriented  area  in  the  reference 
configuration  ©  is  given  by 

dAT  =  eIJKdXJ  dXK .  (83) 


In  the  current  configuration,  this  oriented  area  is  transformed  into  a  new  vector  with  components 
given  by 

dai  =  €ijk  dx2dxk  (84) 


in  s,  or  by 


da-r  =  sj-g  6l^dx]dxk 


(85) 
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in  s.  The  coefficients  of  the  oriented  areas  da*  and  da*  obey  a  simple  tensor  relationship;  they 
represent  the  same  tensor  quantity  in  two  different  coordinate  systems.  The  relationship  between 
these  coefficients  and  the  dA  j  can  be  computed  from  the  definition,  resulting  finally  in 

<9r*  BX1 

da*  =  —  da*  =  \FC  — —  dA^  (86) 

ox1  ox 1 

Volumes  are  computed  in  much  the  same  way.  Reference  volumes  in  &  are  computed  by  the 
Euclidean  scalar  triple  product 

dV  =  dXxUi  •  (dX2U2  x  dX3U3)  =  dX1dX2dX3.  (87) 

Spatial  volumes  in  s  are  computed  likewise,  resulting  in 

dv  =  dx1dx2dx3.  (88) 

The  curvilinear  definition  of  the  vector  product  needs  to  be  invoked  to  compute  volumes  in  s, 
resulting  in 

dfi  =  didui  •  (dT2u2  x  dT3u3)  =  ftg  dT1dT2dT3.  (89) 

These  three  different  expressions  for  two  different  differential  volumes  are  related  by 

dv  =  dv  =  \fC  did  (90) 


since  at  all  times  C  =  g  and  dX 1  =  5~dx\ 

5.3  Three  More  Basic  Relationships 

Before  presenting  the  Maxwell  equations  in  convective  coordinates,  three  more  basic  results  of 
geometry  and  continuum  mechanics  are  needed.  In  the  standard  continuum  mechanics  literature, 
these  relationships  are  often  written  in  reference  coordinates,  but  here  we  translate  them  (where 
needed)  into  convective  coordinates. 


5.3.1  Particle  Velocity  and  Identity  Flux 


The  first  basic  equation  derived  here  relates  the  particle  velocity  components  vl  of  equation  71  to 
the  identity  flux  components  % ;*  of  equation  72  (3).  By  equation  74,  dxl/dt  =  0,  so  that  choosing 
f  =  x1  in  equation  75  gives 


ft 


dxl  ■ 

- 7—V 

dxl 


(91) 
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5.3.2  The  Spatial  Equation  of  Continuity 


The  second  equation  presented  here  is  the  spatial  equation  of  continuity  (1-3,  9).  This 
relationship  is  predicated  on  the  ability  to  compute  the  temporal  derivative  to  the  evolving  metric 
which,  in  turn,  involves  differentiating  the  transformation  tensor  with  respect  to  time.  In 
particular,  the  temporal  derivative  of  the  transformation  tensor  is  given  by 


d  /  dxl  \  dvl  dvl  dxj 
dt  \  dxl  J  did  dxl  Qxl 


(92) 


Given  the  definition  of  the  transformation  determinant,  the  product  rule  of  calculus,  and  the 
vanishing  of  determinants  with  equal  rows,  this  relationship  may  be  used  to  show  that 


d \fg  du*  r: 1 

dt  (lx1  V 9 


(93) 


Given  any  (scalar,  vector,  or  tensor)  quantity  T ,  invoking  the  above  along  with  the  product  rule 
gives  rise  to  the  spatial  equation  of  continuity: 

1  d  (y/gT)  _  87  d(v'T) 
yfg  dt  dt  dxi 


5.3.3  Convected  Time  Derivative 


Finally,  we  have  occasion  to  compute  the  material  derivative  of  a  flux.  Given  a  vector  gl  we  seek 
a  vector  g1  called  the  convected  time  derivative  (3,  9)  of  (f  such  that 

J  gl  dai  =  J  g 1  da*.  (95) 

r  r 


Here,  T  C  M3  is  an  open  surface  in  space;  that  is,  it  is  isomorphic  to  a  finite  part  of  a  plane.  We 
may  compute  the  necessary  quantity  in  several  steps.  First,  by  invoking  equations  30  and  77,  the 
integrand  may  be  brought  into  convective  coordinates  and  the  derivative  moved  inside  the  integral: 


d_ 

dt 


r 


1  dxj 

VF  dl 


€jke  dxkdxe. 


(96) 


Thus,  using  the  product  differentiation  formula  and  equation  94, 


dgl  &  (  k  ,  dxl  d  /  dxl\ 

~dt  +  ~dzf  ^  9  >  +  9  a?  dt  V&9  )  ' 


(97) 
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Using  the  fact  that  f-  is  constant  in  time, 

dxl  d  /  dxl\  dvl 
dxl  dt  \  dxi  )  dxi  ’ 


(98) 


follows  immediately.  Therefore,  inserting  this  relation  into  equation  97,  we  find 


dgl  d(vkgl)  dvl  ■ 

a  =  —  H — - — — - aJ 

■1  EU  1  El  Cl  ..iJ 


dt  dxk  dxi 

which,  upon  relabeling  indices  and  applying  the  product  rule  becomes  the  desired  relationship 


(99) 


dt  dxi 

This  is  the  final  form  of  the  convected  time  derivative. 


(100) 


6.  Maxwell’s  Equations  in  Convective  Coordinates 


6.1  Maxwell’s  Equations  in  Vacuum  in  Spatial  Coordinates 


Before  transforming  the  Maxwell  equations  to  a  convective  frame,  they  must  be  stated  in  the 
usual  spatial  frame.  In  SI  units,  the  basic  quantities  of  electromagnetic  theory  are  total  charge  q 
measured  in  coulombs  per  meter  cubed  (C/m3),  total  current  density  f  measured  in  amperes  per 
meter  squared  (A/m2),  the  electric  field  e,  measured  in  volts  per  meter  (V/m),  and  the  magnetic 
flux  density  If  measured  in  tesla  (T).  We  have  not  demonstrated  that  any  of  these  quantities  is 
tensorial;  therefore,  we  make  these  definitions  at  the  outset  and  do  not  assume  they  have  tensorial 
properties.  In  particular,  we  do  not  change  bases  by  raising  or  lowering  indices.  (These 
definitions  are  the  ones  that  make  sense  later  in  that  they  preserve  the  form  of  the  macroscopic 
Maxwell  equations  in  all  systems.  Therefore,  these  definitions  have  been  made  with  considerable 
hindsight.)  In  vacuum,  the  Maxwell  equations  are 


ijk 


dlf 

dxi 

dek 

dxi 


ijk 


d 


dxi 


j  {Skib1) 


1 

eo  ’ 

0, 

dbl 

~~dt' 

/Ue0  7^  ( ^ej )  +  dof- 


(101) 

(102) 

(103) 

(104) 
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The  quantity  //0  =  4tt  x  10”'  henries  per  meter  (H/m)  is  called  the  permeability  of  free  space ,  and 
the  quantity  e0  =  (//0c2  j” 1  farads  per  meter  (F/m)  is  called  the  permittivity  of  free  space,  where 
c  =  299792458  m/s  is  the  speed  of  light  in  vacuum.  The  electric  field  and  magnetic  flux  density 
described  above  are  basic  in  the  sense  that  the  Lorentz  force  equation 


fi  =  q  (e»  +  €ijkv 3bk) 


(105) 


relates  the  mechanical  force  density  f  (N/m3)  to  these  two  fields,  the  charge  density  q,  and  the 
speed  at  which  the  charge  is  moving  vf  Of  course,  the  current  density  is  related  to  the  charge 
density  by  the  usual  formula  jf  =  qv\  Equation  101  is  known  as  GauB’s  law  for  the  electric 
field;  equation  102  is  known  as  GauB’s  law  for  the  magnetic  field;  equation  103  is  known  as 
Faraday’s  law,  and  equation  104  is  known  as  the  Ampere-Maxwell  law.  Charge  conservation 
follows  immediately  upon  taking  the  divergence  of  equation  104  and  inserting  equation  101, 
resulting  in  the  equation  of  continuity 


djl  dq 
dxl  +  dt 


(106) 


6.2  Maxwell’s  Equations  in  Ponderable  Materials  in  Spatial  Coordinates 


Ponderable,  that  is  massive,  materials  affect  electromagnetic  wave  propagation  because  they  are 
constructed  out  of  charged  particles  like  protons  and  electrons.  Equations  101-104  are  not 
particularly  convenient  in  practice  in  such  media  because  the  charge  and  current  they  reference 
are  the  total  charge  and  current;  that  is,  they  make  no  distinction  between  charges  and  currents 
impressed  externally  to  create  fields  and  those  arising  merely  because  fields  acting  upon  charges 
bound  in  matter  are  acted  upon  by  other  fields.  (Another  way  of  saying  this  is  that  they  are  in 
microscopic  form;  they  consider  all  matter  from  a  corpuscular,  rather  than  continuum,  viewpoint.) 
Bound  charge  densities  arise  either  through  the  creation  of  new  dipoles  or  the  alignment  of 
existing  dipoles  in  response  to  an  external  field.  The  density  of  these  dipoles  can  be  reckoned  as 
a  vector  dipole  moment  per  unit  volume  pl  (C/m2)  related  to  the  bound  charge  density  by 


(107) 


and  this  dipole  density  can  be  related,  in  turn,  to  the  total  field  (13).  Defining  the  free  charge  as 
the  difference  between  the  total  charge  and  the  bound  charge, 


qF  =  q~qB, 


(108) 
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the  GauB  law  for  the  electric  field  becomes 


eo 


_d_ 

dxi 


0 %)  =  »  - 


dp1 
dx i 


This  is  generally  simplified  by  defining  the  electric  displacement  density 


dl  =  e0  +  //, 


(in  C/m2)  and  writing  the  GauB  law  as 


ddl 

d^  =  <lF- 


(109) 


(HO) 


(HI) 


Similarly,  the  total  current  density  may  be  broken  into  its  constitutive  parts: 


f  ^  ?  +  A  +  A 
3  ~Jf+  at  +£ 


TTlk  "T  ^krnnP  V 


(112) 


Here  jp  is  the  free  current  (impressed  and  conducted)  density,  the  term  dp1  /dt  is  the  current 
caused  by  local  increase  or  decrease  of  polarization  charges,  and  the  term  related  to  €kmnPmvn  has 
to  do  with  the  convected  polarization  charges.  Finally,  the  is  the  magnetization  per  unit 
volume  (A/m3),  resulting  from  the  formation  of  magnetic  dipoles  in  reaction  to  the  field  (5,  13). 
Inserting  this  definition  into  equation  104  gives  rise  to  the  macroscopic  Ampere-Maxwell  law 


Mkdh* 

dxi 


Jf  + 


ddl 

W' 


(113) 


where  the  magnetic  field  hk  (in  A/m)  is  defined  as 


hi 


(Sift) 

do 


-  rrii  -  6ijkpPvk. 


(114) 


6.3  Maxwell’s  Macroscopic  Equations  in  Convective  Coordinates 

In  principle,  translating  the  Maxwell  equations  into  convective  coordinate  form  should  be  easy. 
After  all,  if  electromagnetics  is  a  proper  physical  theory,  its  quantities  should  change  coordinate 
systems  as  tensors,  and  formulas  can  be  written  immediately.  Unfortunately,  that  basic  physical 
observations  prohibit  this  is  clear  upon  momentary  reflection:  An  observer  holding  a  charge  sees 
an  electric  field,  but  another  observer  moving  with  respect  to  the  first  sees  a  current  and  therefore 
a  magnetic  field.  This  implies  immediately  that  the  electric  and  magnetic  fields  must  be  part  of 
the  same  physical  entity,  as  indeed  they  are  in  the  theory  of  relativity  (10,  11,  13).  Despite  this, 
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nonunique  theories  may  be  constructed  via  other  means;  this  is  the  approach  taken  here  to 
illustrate  the  convective  form  of  reference  3.  This  work  aimed  to  preserve  the  macroscopic 
Maxwell  equations.  Though  other  approaches  are  possible,  we  do  not  consider  them  here. 


To  convert  equation  113  to  convective  coordinates,  we  integrate  both  sides  over  an  arbitrary 
surface  and  apply  Stokes’s  theorem  to  find 


J  =  (f  didx1. 

r  av 


(115) 


Here,  the  covariant  components  of  the  spatial  directed  differential  area  are  given  by  equation  84, 
and  the  right-hand  side  indicates  a  line  integral  around  the  closed  path  <9T  surrounding  T  in  spatial 
coordinates.  (The  path  c)T  has  direction  related  to  that  of  the  normal  to  T  by  the  right-hand  rule. 
The  circle  notation  on  the  path  integral  indicates  that  the  path  is  necessarily  closed.)  We  can 
simplify  the  right-hand  side  of  equation  113  using  the  convected  time  derivative  defined  in 
equation  100.  Specifically,  using  equation  100  to  substitute  for  the  partial  derivative  of  the 
electric  displacement,  and  recognizing  the  meaning  of  that  derivative,  equation  113  can  be  written 
in  the  form 

y  (hi  -  eijkvJ dk)  dxl  =  ^  J  d1  dat  +  J  (^f  ~  v%^^j  da*.  (116) 

dr  r  r 

Because  the  time  derivative  in  this  equation  is  the  material  time  derivative,  the  integrals  can  be 
converted  to  convective  coordinates  using  the  techniques  from  the  previous  sections. 

Recognizing  the  free  charge,  this  procedure  results  in  the  final  integral  form  of  the 
Ampere-Maxwell  law  in  convective  coordinates, 

/f)  r%  _  A  /•  r  .  0;p 

0 i  -  eijkV3dk)  -^dxl  =  —  J  d'—ddt  +  J  (j'F  -  qFvl)  —  ddh  (117) 

sr  r  r 


and  Stokes’s  theorem  gives  its  final  differential  form: 


d_ 


{hi  -  ( ijkiJ (I1') 


dxl 
dx * 


d  /  ^  d:F 
dt  \  dx{ 


+  {jf  ~  Qfv%) 


dx1 

dxi 


(118) 


Because  we  want  the  form  of  this  equation  preserved  under  the  coordinate  change,  we  insist  that 
this  equation  be  interpreted  as  the  convective  Ampere-Maxwell  law 


1 

7F 


ddi 

dt 


+  ic, 


(119) 
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so  that  we  define  the  convective  magnetic  field  by 


hi  ( hi  ( ]  j\;  (1  ^ 

the  convective  displacement  density  by 


dxl 
<957  ’ 


dl  =  —dl 
dx'  ’ 


and  the  conducted  convected  current  by 


io  =  (A  -  ifV)  Q-. 


(120) 


(121) 


(122) 


(This  last  is  called  conducted  current  and  not  free  current  because  the  only  charge  moving  relative 
to  the  convected  frame  must  be  that  conducted  through  it.)  Faraday’s  law  is  converted  in  the 
same  way  and  becomes 


1  mdh  =  _<& 

yfg  dxJ  d  t  ’ 


(123) 


where 


and 


e~i  =  {et  +  eijkvjbk )  — , 


(124) 

(125) 


The  divergence  equations  can  be  converted  even  more  simply.  Integrating  both  sides  of 
equation  111  over  an  arbitrary  volume  and  applying  the  divergence  theorem  immediately  gives 
the  convective,  integral  form  of  the  GauB  law  for  the  electric  field: 


/ 


•  dr1 

dl—dd~t  = 
dxl 


(126) 


Here,  Q  C  M3  is  a  volume  and  0i l  is  the  surface  enclosing  it.  The  circle  around  the  integral 
symbol  used  on  the  left-hand  side  of  this  equation  implies  a  surface  that  is  necessarily  closed.  If 
we  define 

qF  =  qF,  (127) 


and  apply  equation  121  and  the  divergence  theorem  again,  the  convective  differential  form  of  the 
electric  GauB  law  obtains: 


(128) 
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The  magnetic  GauB  law 


VF5ii(V®6‘)  °’ 


(129) 


is  derived  similarly. 

6.4  Maxwell’s  Microscopic  Equations  in  Convective  Coordinates 

Many  of  the  same  techniques  can  be  applied  to  transform  the  microscopic  Maxwell  equations  to 
convective  coordinates.  That  this  needs  to  be  discussed  at  all  springs  from  the  fact  that  these 
equations  directly  incorporate  the  contravariant  components  of  the  electric  field  vector  and  the 
covariant  components  of  the  magnetic  flux  density,  in  spite  of  the  fact  that  these  cannot  be  derived 
from  the  usual  process  of  raising  or  lowering  indices.  (That  this  is  the  case  is  already  clear  from, 
for  instance,  equation  124;  the  coordinate  change  formula  for  the  electric  field  should  not  involve 
the  magnetic  field  if  the  electric  field  is  a  true  independent  tensor  quantity.) 

On  the  other  hand,  two  of  the  microscopic  Maxwell  equations  are  identical  to  their  macroscopic 
counterparts  because  they  are  homogeneous.  These  two  equations,  the  GauB’s  law  for  the 
magnetic  field  (equation  129)  and  the  Faraday  law  (equation  123)  therefore  are  not  altered.  The 
two  equations  containing  sources,  however,  are  more  complicated  and  we  turn  to  them  now. 

6.4.1  GauB’s  Law  for  the  Electric  Field 


Integrating  GauB’s  law  for  the  electric  field  in  microscopic  form  (equation  101)  over  an  arbitrary 
volume,  and  using  Stokes’s  theorem  to  convert  the  field  integral  to  a  surface  integral  yields 


eo  <P  S^e.j  da* 


on 


(130) 


Using  equations  86  and  90,  these  areas  and  volumes  can  be  converted  to  convective  coordinates 
resulting  in 


eo 


Sij— _ e  da~ 

dx'  31 


q  dv. 


(131) 


on  ci 

Using  the  divergence  theorem  again  to  convert  the  left-hand  side  of  this  equation  back  into  a 
volume  integral,  and  defining  the  convected  total  charge  in  the  obvious  way  (i.e.,  q  =  q ),  we  find 


r  i  d 
J  y/g'dx1 


■  .  c)t^ 
dxi  3 


di) 


ci 


(132) 
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Finally,  recognizing  the  arbitrariness  of  the  integration  volume,  and  substituting  the  expression 
for  the  electric  field  from  the  inverse  of  equation  124  (while  invoking  the  velocity  transformation 
of  equation  91),  we  find  the  final  microscopic  form  of  the  GauB  law  for  the  electric  field: 


eo 


1  d 

\fg  dxl 


=  q- 


(133) 


We  can  derive  the  definition  of  the  electric  displacement  density  in  terms  of  the  electric  field  and 
polarization  density  from  this  equation.  By  expressing  the  total  charge  in  terms  of  its  free  and 
bound  components  in  equation  131  (while  substituting  the  expression  for  electric  field  from 
equation  133),  then  proceeding  on  the  left-hand  side  as  before,  we  can  write 


eo 


1  d 
\fg  dxl 


\fa 


(ej  +  y/ff' Wkrh$k~b 


dv=  Udv-  j8£iv. 


(134) 


Now,  following  the  same  steps  as  before  to  convert  the  rightmost  integral  first  to  a  surface  integral 
and  back,  changing  to  the  convected  coordinate  system  along  the  way,  the  right-hand  side 
becomes 

/  «F  d"  -  /  t?dF  =  /  «F  ds  “  /  f  (^p‘) ds-  <135) 

n  no 

Thus,  defining  the  convected  polarization  density  through  the  standard  tensorial  transformation 


.  dxl  i 

V  =  -7T-V 

ox1 


(136) 


we  may  write  that  in  the  convected  coordinates 

_  _  Id 


q  =  qF 


yjg1  dxl 


(137) 


so  that 
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1  d 
\fg  dxl 
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r  (e, + VF + 
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du  =  J  qF  du. 
n 


(138) 


Comparison  of  this  equation  with  equation  128,  coupled  with  the  realization  that  the  volume  of 
integration  is  arbitrary,  leads  to  the  final  form  of  the  convected  constitutive  law: 


d™  =  eo~g™  (e,+  y/?> + 


■ P " 


(139) 


This  is  not  of  the  form  d  =  e0e  +  P  seen  in  s.  Therefore,  if  we  are  to  believe  this  equation,  we 
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must  believe  that  the  laws  of  physics  depend  on  the  observer,  in  direct  contradiction  to  the 
principle  of  relativity.  The  reader  may  object  that  this  particular  equation  results  only  from  our 
decision  to  force  the  macroscopic  Maxwell  equations  to  have  the  same  form  in  different 
coordinate  systems.  This,  of  course,  is  true,  but  if  we  attempted  to  force  the  constitutive  relation 
to  hold  in  all  systems,  we  would  find  the  macroscopic  Maxwell  equations  would  be  of  different 
form  in  each  system.  We  conclude  from  this  that  electromagnetics  is  incompatible  with  the 
Galilean  relativity  assumed  in  standard  Newtonian  physics. 

6.4.2  The  Ampere-Maxwell  Law 

We  finally  turn  to  the  convective  form  of  equation  104,  the  microscopic  Ampere-Maxwell  law. 
Integrating  this  equation,  and  using  Stokes’s  theorem  and  the  definition  of  the  convected  time 
derivative  (equation  100),  we  find 

/  ib  -  ^  <^1  d*‘ = 4  /  **  ** + /  V  -  4?  "1  d“- <140) 

dr  r  r 

The  integrals  are  now  converted  over  to  convective  coordinates  using  the  usual  formulas,  giving 

/I  1  f)ri 

-bi-e0eijkvj(Skmem)  dx* 

or 

ci  f  f  a  a^a 

=  £»d? ./  dSi  + ./  Y  ~  £“a?  ^  1  a?  da"  <141) 

r  r 

From  the  last  term  of  this  equation,  we  immediately  recognize  the  convective  total  current 

arl 

f  =  a?  O'  - 1«‘)  ■  a«> 

Now,  the  left-hand  side  of  equation  141  can  be  manipulated  algebraically  to  show  that 

— bi  -  e0 eijkv3  ( 5kmem )  =  —g~tmbm  +  e0 \fg 'e~^3gkrn  (e^  +  a Jq .  (143) 

ji0  J  /^o  v  ' 

Inserting  these  relations  into  equation  141,  and  once  again  invoking  the  inverse  of  equation  124 
and  recognizing  the  arbitrariness  of  the  integration  volume,  gives  the  final  form  of  the 


29 


microscopic  Ampere-Maxwell  law: 


^h:!>’  +  '“A'*",,,’?”')'1'  (4  -  | 

=  (?J  +  ^ W^)]  +  f'  <144> 

This  result  is  inconsistent  with  the  principle  of  relativity;  it  looks  nothing  like  its  spatial  version. 

Finally,  we  need  to  derive  the  magnetic  constitutive  relations  in  convective  coordinates  from  this 
equation.  We  proceed  in  a  manner  similar  to  that  for  the  electric  relationship.  To  begin,  we  can 
separate  the  current  and  charge  into  their  component  parts  using  equations  112  and  108,  and 
insert  these  definitions  into  the  integral  of  equation  144.  Thus,  equation  144  is  written  as 


The  conduction  current  is  immediately  recognizable  here,  and  so  can  be  simplified  out  of  the 
messy  last  integral.  The  resulting  expression  reads 


/  +  e°^ eifkft9ki  (e/  +  \[g  dx ' 

or 

=  e°^i  f  ^  (gj  +  ^ ejfsVbs^  da* 

f 

+  /  I?  {  fr  +  ^  mk’j  +  (€kimPtvm)  +  pdyi]  d5*  +  J ?cdah  (146) 
r  r 

where 

/o  =  Q~fc  =  Of  ~  9fV1)  =  Jp  +  qFft.  (147) 

Next,  by  transforming  the  Levi-Civita  symbol  and  invoking  Stokes’s  theorem,  we  can  show  that 


r 


dx1  ijkd mk 
dx 1  dxi 


er 


(148) 
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from  which  we  gather  that 


dxk 

dxk 


mk. 


(149) 


Finally,  we  invoke  the  definition  of  the  convected  time  derivative  (equation  100)  to  simplify  the 
terms  involving  the  polarization: 


dxl 

dxi 


dt  dxi 


f)ry) 

(ekempvm)  +  —v 


dxi 


da,- 


(150) 


Substituting  all  of  these  equations  back  into  equation  146  finally  yields 


ar 


[— a4>1  +  eo VF '^fki>3gu  +  \fg  -  m 


dxl 


=  eo 


dt 


f/'3  (<%  +  eFStjjrbs)  +pl  da-t+  fcda^.  (151) 


Recognizing  the  expression  for  dl  buried  in  here,  we  find  that  the  constitutive  relationship 

h  =  —9ub+e o  e-lfk^gM  (eE  +  y/F -  m  (152) 

IM)  v  / 

recovers  equation  1 19.  Notice  that  once  again,  this  equation  does  not  resemble  its  spatial  form. 
Finally,  a  careful  accounting  of  the  dispensation  of  all  of  the  terms  in  the  above  demonstrates  that 


A  _  Jfkdfnk  d $ 
3  ~Jce  ~W  +  ~dt 


(153) 


The  convective  contribution  of  the  polarization  density  present  in  the  spatial  domain  is  absent 
here  as  it  is  bound  to  the  material  that  is  not  moving  in  the  this  frame. 


7.  Conclusions 


This  report  has  documented  the  formulation  of  physical  laws  in  moving  reference  frames,  with 
extra  attention  paid  to  the  Maxwell  equations  of  electromagnetic  theory.  The  first  part  of  the 
report,  making  up  its  bulk,  describes  how  the  expression  of  physical  theories  change  with 
changing  coordinate  systems.  In  particular,  this  portion  of  the  work  discusses  what  sort  of 
relationships  physical  quantities  must  obey  so  that  the  resulting  theory  is  consistent  when  the 
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coordinate  system  that  forms  the  basis  for  its  description  is  curvilinear,  nonorthogonal,  and 
time-varying.  In  particular,  the  discussion  shows  that  two  different  sets  of  basis  vectors  are 
convenient  for  the  expression  of  physical  law  in  nonorthogonal  systems:  A  co variant  set  (called 
u,  in  the  text),  in  which  each  basis  vector  points  in  the  direction  of  variation  of  a  given  coordinate 
holding  all  other  coordinates  fixed,  and  a  contravariant  set  (called  fi'  in  the  text)  pointing 
orthogonal  to  the  constant  coordinate  surfaces.  In  orthonormal  systems,  these  bases  are  identical, 
but  in  nonorthogonal  systems  they  are  merely  biorthogonal. 

These  different  bases  are  chosen  merely  because  they  are  useful  for  the  expression  of  physical 
laws,  and,  in  general,  have  no  deep  physical  significance  in  and  of  themselves.  Changing  from 
one  system  to  the  other  is  merely  a  mathematical  change  of  basis;  for  the  underlying  physical 
theory  to  be  meaningful,  such  a  transformation  must  be  physically  irrelevant.  (A  theory  that  said 
otherwise  would  imply  that  the  world  is  affected  by  our  description  of  it!)  When  the  Maxwell 
equations  are  formulated  in  the  material  coordinate  systems  of  continuum  mechanics  this  simple 
consistency  requirement  fails,  and  for  this  reason,  the  literature  is  filled  with  a  profusion  of 
contradictory  theories  about  how  the  Maxwell  equations  should  be  expressed  in  the  coordinate 
systems  used  by  continuum  mechanics  other  than  the  spatial. 

The  second  “half”  of  this  report  seeks  to  clarify  the  expression  of  the  Maxwell  equations  in 
coordinate  systems  relevant  to  continuum  mechanics.  In  so  doing,  it  subtly  illustrates  two 
primary  points: 

•  The  consistency  or  inconsistency  of  a  theory  cannot  be  determined  by  its  expression  in 
material  or  “Lagrangian  ”  coordinates,  but  must  be  determined  in  a  coordinate  system 
describing  the  same  physics  as  the  spatial  system,  such  as  the  convective  coordinate  system. 

•  Under  standard  Newtonian  physics,  invoking  the  Galilean  transformation  for  moving 
systems,  the  convective  form  of  the  Maxwell  laws  cannot  be  consistently  formulated, 
resulting  in  the  confusion  in  the  literature.  Consistent  formulation  requires  the  theory  of 
relativity  and  is  the  subject  of  the  sequel  to  this  report. 

We  discuss  each  of  these  points  in  turn. 

The  first  point  is  perhaps  most  clearly  understood  as  it  results  from  the  definition  of  material 
coordinates  as  merely  labels  for  material  points;  they  function  as  a  way  of  tracking  individual 
particles  on  their  journey  through  space  and  time.  For  this  reason,  not  only  need  they  not 
correspond  to  the  state  of  the  body  in  question  at  the  initial  time,  they  need  not  match  the  physical 
arrangement  of  the  body  at  any  time  whatsoever.  They  need  merely  represent  a  potential  state  of 
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the  body  so  that  the  mapping  describing  its  motion  retains  some  simple  topological  consistency. 
As  a  mere  labeling  of  points,  there  is  no  need  for  any  physical  theory  to  be  satisfied  in  material 
coordinates;  the  configuration  they  describe  is  a  fiction.  This  leaves  the  theorist  at  a  loss  in 
distinguishing  between  different  expressions  of  the  theory  in  the  material  frame,  since  there  is  no 
reason  for  preferring  one  transformation  to  any  other. 

In  this  regard,  we  also  note  that  for  precisely  the  same  reason,  the  experimentalist  cannot  resolve 
the  issue.  The  “value  of  the  electric  field  in  material  coordinates”  is  not  something  that  can  be 
measured  since  it  refers  to  a  coordinate  system  that  may  never  have  existed  in  the  lab.  This  is  not 
a  mere  matter  of  accessibility  either;  there  is  simply  no  field  that  can  be  measured  or  even  suitably 
manipulated  a  posteriori  into  a  quantity  that  can  be  coherently  called  the  material  electric  field. 

The  second  point  is  made  clear  from  section  6.  By  forcing  the  macroscopic  Maxwell  equations  to 
have  the  same  form  in  convective  coordinates,  some  variables  are  forced  to  be  covariant  and 
others  are  forced  to  be  contravariant.  While  the  Maxwell  equations  still  stand  in  a  recognizable 
form  with  the  proper  choices,  the  resulting  constitutive  laws  are  bizarre.  The  equations 
transforming  fields  between  frames  are  coupled,  which  should  not  be  the  case  if  the  electric  and 
magnetic  fields  are  independent  physical  entities.  Worse  still,  indices  cannot  be  “raised”  or 
“lowered”  in  the  usual  way,  a  situation  literally  tantamount  to  saying  that  the  physics  of  the 
situation  depends  on  a  choice  of  basis,  i.e.,  on  an  arbitrary  issue  of  mathematical  depiction. 

These  problems  (with  the  description  in  convective  coordinates)  spring  from  the  constancy  of  the 
speed  of  light  relative  to  all  observers  predicted  by  the  Maxwell  relations  but  not  by  Newtonian 
physics,  and  they  cannot  be  remedied  without  appeal  to  the  special  theory  of  relativity. 
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